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Abstract 

Assume that Au — f, (1) is a solvable linear equation in a Hilbert space H, A 
is a linear, closed, densely defined, unbounded operator in H, which is not boundedly 
invertible, so problem (1) is ill-posed. It is proved that the closure of the operator 
(A* A + aI)~ 1 A*, with the domain D(A*), where a > is a constant, is a linear bounded 
everywhere defined operator with norm < 1. This result is applied to the variational 
problem F(u) := \\Au — f\\ 2 + a||w|| 2 = min, where / is an arbitrary element of H, 
not necessarily belonging to the range of A. Variational regularization of problem (1) is 
constructed, and a discrepancy principle is proved. 
1. Introduction 

The main results of this paper are formulated as Theorems 1 and 2 and proved in 
Sections 1 and 3, respectively. In Section 1 we formulate Theorem 1 which deals with 
a linear, unbounded, closed, densely defined operator A. In Section 2 this operator is 
assumed not boundedly invertible and the problems arising in the study of variational 
regularization of the solution to the equation 

\ Au = f, (1) 

are studied, where A : H — > H is a linear, unbounded, closed, densely defined, not 
boundedly invertible operator on a Hilbert space H with domain D(A) and range R(A). 
Since A is densely defined and closed, its adjoint A* is a closed, densely defined linear 
operator. The operators T = A* A and Q = A A* are nonnegative, selfadjoint, densely 
defined in H operators (see [1]), the operator T a := T + al, (I is the identity operator 
and a > is a constant) is boundedly invertible, i.e., its inverse is a bounded linear 
operator, defined on all of H, with norm < i. It is easy to check that the operator 
A*Q~ 1 is bounded, defined on all of H, and H^Q" 1 !! < We assume in Section 2 

that the operator A is not boundedly invertible, in which case problem (1) is ill-posed. 
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We are interested in the operator S := S a := T~ 1 A* defined on a dense set D(A*). 
The reasons for our interest will be explained soon. The product of an unbounded closed 
operator (A* in our case) and a bounded operator {T~ l in our case) is not necessarily 
closed, in general, as a simple example shows: if A = A* > is unbounded selfadjoint 
(and consequently closed) operator and B = (I + A)' 1 is a bounded operator, then the 
operator BA with the domain D(A) is not closed. Its closure is a bounded operator 
defined on all of H. This closure is uniquely defined by continuity. 

Lemma 1. If A is a linear, closed, densely defined, unbounded operator in H, and 
B is a bounded linear operator such that R(B*) C D(A), then the operator BA* with the 
domain D(A*) is closable. 

Proof. To prove the closability of BA* one has to prove that if u n — > and BA*u n — > 
w, then w = 0. Let h be arbitrary. Then B*h belongs to D(A). Therefore 

(w, h) = \im(BA*u n , h) = \im(u n , AB*h) = 0. 

Thus, w = 0. □ 
In our case B = T~ x = B* and R(T~ l ) C D{A). By Lemma 1, the operator T~ l A* is 
closable. The operator AT~ l is bounded, defined on all of H, with norm < Indeed, 

by the polar decomposition one has A = UT 1 / 2 , where U is an isometry, so \\U\\ < I. 
Thus, \\AT^\\ < IITVX- 1 !! = sup s > gS = ^. 

Lemma 2. The operator S := T^A* with domain D(A*) has the closure S = S** , 
which is a bounded operator defined on all of H, with norm < an d S* = AT^ 1 is a 
bounded operator defined on all of H. 

Proof. The operator S is densely defined. By Lemma 1 it is closable, so the operator 
S* is densely defined. Let us prove that S* = AT' 1 . Let h e H be arbitrary. We have 

(T- l A*u,h) = (A*u,T^h) = (u,AT^h). 

This implies that D(S*) = H and (T^A*)* = AT' 1 . We have used the relation 
RiT^ 1 ) C D{A). Let us check this relation. Let g G R{T~ l ), then g = T~ l h and 
h = T a g. Thus, g G D(T) C D(A), as claimed. Lemma 2 is proved. □ 
From Lemmas 1 and 2 we obtain the following result. 

Theorem 1. Let A be a linear, closed, densely defined, unbounded operator in a 
Hilbert space. Then the operator S = T^A* with domain D(A*) admits a unique closed 
extension S defined on all of H , with the norm < 

Why should one be interested in the above theorem? 

The answer is: because of its crucial role in the study of equation (1) and of variational 
regularization for equation (1). The corresponding theory is developed in Section 2 and 
proofs are given in Section 3. 

2. Variational regularization 

Assumption A: We assume throughout that A is linear, unbounded, densely defined 
operator in H, and that A is not boundedly invertible, so problem (1) is ill-posed. We 
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assume that equation (1) is solvable, possibly nonuniquely, that f ^ 0, and denote by y 
its unique minimal-norm solution, yA-N , where N = N(A). 

This assumption is not repeated below, but is a standing one throughout the rest of 
this paper. 

Assume that \\fg — f\\ < 5, where fg is the "noisy" data, which are known for some 
given small 5 > 0, while the exact data / are unknown. The problem is to construct a 
stable approximation ug to y, given the data {A, S, fg}. Stable approximation means that 
lim^o IK _ y\ \ = 0. Variational regularization is one of the methods for constructing 
such an approximation. 

If A is bounded, this method consists of solving the minimization problem 

F(u) := F Qj5 = \\Au- f 5 \\ 2 + a\\u\\ 2 = min, (2) 

and choosing the regularization parameter a = a(5) so that lim^o^ = y, where us := 
u a (s),s- It is well known and easy to prove that if A is bounded, then problem (2) has 
a unique solution, u a> g = T~ 1 A* fg, which is a unique global minimizer of the quadratic 
functional (2), and this minimizer solves the equation T a u a ^ = A* fg. The last equation 
does not make sense, in general, if A is unbounded, because fg may not belong to D(A*). 
This is the difficulty arising in the case of unbounded A. In this case it is not a priori clear 
if the global minimizer of functional (2) exists. We prove that this minimizer exists for 
any fg G H, that it is unique, and that there is a function a = a(S) > 0, lim^o a ($) — 0, 
such that lim^o u a(S),s = y, so the element ug := u a (s),s is a stable approximation of 
the unique minimal- norm solution to (1). Theorem 1 allows one to define the element 
T~ l A*fg for any fg, and not only for those fg which belong to D(A*). We also prove for 
unbounded A a discrepancy principle in the following form. Let ug >a solve (2). Consider 
the equation for finding a = a(S): 

\\Au a , s -f s \\ = C6, C>1, \\f s \\>C6, (3) 

where C is a constant. Equation (3) is the discrepancy principle. We prove that equation 
(3) determines a(S) uniquely, a(5) — > as S — > 0, and ug := us >a (S) ~ ¥ V as S — > 0. This 
justifies the discrepancy principle for choosing the regularization parameter (see [2] for 
various forms of the discrepancy principle). 
Let us formulate the results. 

Theorem 2. For any f e H the functional F{u) = \ \Au — f\ | 2 + a\ \u\ \ 2 has a unique 
global minimizer u a = A*Q~ 1 f, where Q = AA* , Q a := Q + al , a > is a constant, 
and 

A*Q~ l f = T~'A*f, (4) 
where T~ 1 A* is the closure of the operator T~ l A* defined on D(A*). If f e R{A), then 

Urn JK- /|| = 0, (5) 

where u a is the unique global minimizer of F{u) and y is the minimal-norm solution to 
(1). If \\fs — /|| < 5 and u a< g = A*Q~ 1 fg, then there exists an a(5) > such that 

lim|K-/|| = 0, lima(<5) = 0, ug := u a(5)j5 . (6) 

8— >0 <5— >0 
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Equation (3) is uniquely solvable for a, and for its solution a(S) equation (6) holds. 
In Section 3 proof of Theorem 2 is given. 
3. Proof of Theorem 2 

3.1. For any h G D(A) let u a := A*Q~ l . One has 

F(u a + h) = F(u a ) + \\Ah\\ 2 + a\\h\\ 2 + 2^[(Au a - f,Ah) + a(u a , h)}, (7) 

and 

{Au a - f, Ah) + a(u a , h) = (QQ' 1 f - f, Ah) + a(u a , h) = -a(Q" 1 /, Ah) + a(u a , h) = 

-a[(A*Q- 1 f + u a ,h) = 0. (8) 

From (7) and (8) it follows that u a is the unique global minimizer of F(u). 

3.2. Let us prove (4). If (4) holds on a dense in H linear subset D(A*), then it holds 
on all of H by continuity because A*Q~ 1 is a bounded linear operator, defined on all 
of H, with norm < so that the closure of the operator T~ l A* defined on D(A*), 

is a bounded operator, defined on all of H, with norm < Indeed, let / G D(A*), 

g := Qa l .fi so Qa9 = f and g G D(A*AA*). Therefore equation (4) is equivalent 
to A*Q a g = T a A*g, or A*AA*g + aA*g = A*AA*g + aA*g, which is an identity. If 
/ G D(A*) (so that g G D(A*AA*)), then the above formulas are justified and one can 
go back from the identity A*AA*g+aA*g = A*AA*g+aA*g, valid for any g G D(A*AA*), 
define / = Q a g, (this / belongs to D(A*) because Qg G D(A*)), and get (4). 
Note that if A were bounded, then one would have the identity 

A*(f)(Q) = (p(T)A*, T = A* A, Q = AA* , (9) 

valid for any continuous function <p. Indeed, if is a polynomial, then (9) is obvious 
(for example, if <f>{s) = s, then (9) becomes A*(AA*) = (A*A)A*). If is a continuous 
function on the interval [0, ||^4|| 2 ], then it is a limit (in the sup-norm on this bounded 
interval) of a sequence of polynomials (Weierstrass' theorem) , so (9) holds. In our problem 
A is unbounded, so are Q and T, and (f>{s) = (with a = const > 0) is a continuous 
function on an infinite interval [0, oo). Linear unbounded operators do not form an 
algebra, in general, because of the difficulties with domain of definition of the product of 
two unbounded operators (the product may have the trivial domain {0}). That is why 
formula (4), which is a particular case of (9) for bounded operators, has to be proved 
independently of this formula. 

3.3 Let us prove (5). If / G R(A), then / = Ay, where y~LN is the minimal-norm 
solution to (1). We have u a — y = T~ l Ty — y = —aT~ x y and 

f°° a 2 

lim n I \ aT a x y\ I 2 = lim n / 7 — —^ d ( E sV, v) = I \Pnv\ I 2 = o, 

where E s is the resolution of the identity of the selfadjoint operator T and Pjy is the 
orthogonal projector onto iV = N(A), so ||-Pjvy|| = because y-LN. 
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3.4. Let us prove (6). We have 

I \us - y\ I < I \ u s -u a \\ + \ \u a - y\ | := h + I 2 - 
We have already proved that lim^o I 2 = 0, because lim^o a(S) = 0. Let us estimate I±: 

5 



u s - u c 



A*Q-J 5) (fs-f)\\< 
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Thus, if lim^o a(5) = and lim^o 27s = ^' ^ nen (^) holds. 
3.5. Finally, let us prove 
The discrepancy principle: 

Equation (3) is uniquely solvable for a and its solution a (5) satisfies (6). 
The proof follows the one in [2], p. 22. One has 

g(a,5) := \\Au a , s - fs\\ 2 = \\QQ?f, - fs\\ 2 = a 2 f°° d f s [ sJs J = C 2 5 2 , (10) 

Jo \ s + a ) 

where S s is the resolution of the identity of the selfadjoint operator Q. The function 
g(a) := g(a, 5) for a fixed 5 > is continuous, strictly increasing on [0, 00) and g(oo) > 
C 2 S 2 while g(0) < 5 2 , as we will prove below. Thus, there exists a unique ct(8) > 0, 
such that g(a(5), S) = C 2 S 2 , and lim^o a ($) — because g(a, 5) > for a 7^ and any 
5 G [0, 5 ), provided that ||/|| 7^ 0, which we assume. Here 6 > is a sufficiently small 
number. 

Let us prove the two inequalities: g(oo) > C 2 S 2 and g(0) < S 2 . We have 

PCX) 

g(oo) = / d(£ s f s ,fs) = \\fs\\ 2 >C 2 5 2 , 
Jo 

because of the assumption \ \fs\ \ > CS. Also 

g(0) = \\PN(Q)f8\\ 2 <S 2 . 

Indeed, f s = f + (f 5 - /). The element / G R(A), so f±N(A*) = N(Q). Therefore 
Pn(Q)Js = PN(Q)(fs ~ /)■ Consequently, 

\\PN { Q)fs\\ < \\PN(Q)(fs ~ f)\\ < \\fs-f\\<S. 

Let us prove the limiting relation lim^o \ \ u s ~ f\ \ — 0. We have 

F a{ s),sM = \\Au s - fs\\ 2 + a(5)\\u s \\ 2 < F am (y) < 5 2 + a(S)\\y\\ 2 . (11) 

Since ||Au,5 — fs\\ 2 = C 2 b 2 > 5 2 we conclude from (11) that \ \us\ \ < \\y\ \ for all 5 G [0, 5q). 
Thus we may assume that u$ — 1 z as 5 — > 0, where — ^ denotes weak convergence in H. 
Since lim^ fs = f, we conclude from \ \Aus — fs\\ = CS that lim^o \ \Aus — f\\ — 0. 
This implies Az = f. Indeed, for any h G D(A*) one has 

(f,h) = \im(Au s ,h) = (z,A*h). 

8— »0 



Therefore Az = f. Since \\us\\ < ||y||, we have lim^ol \ u s\ I < \\v\\- From u$ — z we 
obtain ||z|| < lim g _, ||^|| < Thus, \\z\ \ < \\y\\. Since the minimal- norm solution to 
(1) is unique, it follows that z = y. Thus, u$ — 1 y and lim^ol |m<s| | < ||y|| < lim^ | \us\ |. 
This implies lim^o H^H = Consequently, lim^o — y\ \ — 0, because \ \us — y\\ 2 = 
I Kl ? + I |y| I 2 - 23?(u 5 , y) -»• as 5 -> 0. 

Theorem 2 is proved. □ 
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